source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/BGH/BGH Annual Gini.R')
g = aggregate(d, list(d$Year), FUN = ineq)
View(g)
boxplot(g$Pay)
View(g)
boxplot(g, ylim = c(0.1,0.3))
boxplot(g$Pay, ylim = c(0.1,0.3))
boxplot(g$Pay, ylim = c(0.15,0.25))
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
setwd("~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/BGH")
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
setwd("~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/BGH")
library(readr)
library(ineq)
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
test = subset(d, d$Year = 85)
test = subset(d, d$Year == 85)
test = subset(test, test$Level == 1)
hist(log(test))
hist(log(test$Pay))
library(readr)
library(ineq)
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
d$Year = paste("19", d$Year, sep = "")
d$Year = as.numeric(d$Year)
y = sort(unique(d$Year))
l = length(y)
final = NULL
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
d = read_csv("BGH.csv")
d = read_csv("BGH.csv")
d = read.csv("BGH.csv")
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
d$Year = paste("19", d$Year, sep = "")
d$Year = as.numeric(d$Year)
View(d)
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
View(d)
d = read_csv("BGH.csv")
d = data.frame(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
d$Year = paste("19", d$Year, sep = "")
d$Year = as.numeric(d$Year)
library(data.table)
d = data.table(Year = d$year, Level = d$level, Pay = d$salary)
d = read_csv("BGH.csv")
d = data.table(Year = d$year, Level = d$level, Pay = d$salary)
d = read_csv("BGH.csv")
d = data.table(Year = d$year, Level = d$level, Pay = d$salary)
d = na.omit(d)
d$Year = paste("19", d$Year, sep = "")
d$Year = as.numeric(d$Year)
View(final)
View(final)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
summary(d$Pay)
plot(d$Year, d$Pay, log = "y")
plot(d$Level, d$Pay, log = "y")
plot(d$Level, d$Pay)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
plot(d$Level, d$Pay)
plot(d$Level, d$Pay, log = "y")
View(final)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/BGH Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
test = d[,list(
Year = head(Year, 1),
Level = head(Level,1)
),
by = .(Year, Level)       ]
View(test)
View(test)
test = d[,list(
N = length(Pay),
Mean = mean(Pay)
),
by = .(Year, Level)       ]
View(test)
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = ineq(Pay)
),
by = .(Year, Level)       ]
View(test)
list(
N = length(Pay),
Mean = mean(Pay),
Gini = ineq(Pay)*N
)
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = ineq(Pay)*length(Pay)/(length(Pay) - 1)
),
by = .(Year, Level)       ]
View(test)
?ifelse
gini.adjust = function(x){
ifelse(length(x) == 1, 0, ineq(x)*length(x)/(length(x) -1))
}
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = gini.adjust
),
by = .(Year, Level)       ]
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = gini.adjust(Pay)
),
by = .(Year, Level)       ]
View(test)
mean(test$Gini)
mean(final$Gini)
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = gini.adjust(Pay)
Source = "BGH"
),
by = .(Year, Level)       ]
test = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = gini.adjust(Pay),
Source = "BGH"
),
by = .(Year, Level)       ]
View(test)
plot(final$Gini, test$Gini)
result = d[,list(
N = length(Pay),
Mean = mean(Pay),
Gini = gini.adjust(Pay),
Source = "BGH"
),
by = .(Year, Level)       ]
result = result[order(Year, Level), ]
View(result)
plot(final$Gini, result$Gini)
sd(1)
sd(5)
plot(result$Mean, final$Mean)
d = read.csv("BGH.csv")
d = data.table(Year = d$year, Level = d$level, Pay = d$salary) %>% na.omit()
d$Year = paste("19", d$Year, sep = "")
d$Year = as.numeric(d$Year)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Level Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Promotions Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/BGH/Promotions Analysis.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Case Studies/Dohmen/Data/Dohmen Firm Gini.R')
